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o ■ 

A compact binary in an eccentric orbit radiates gravitational waves (GWs) at all in- 
. teger harmonics of its orbital frequency. In this study, we investigate the effect of orbital 

D ' eccentricity on the expected gravitational background radiation (GWBR) from supermassive 

PL| ; black hole (SMBH) binaries in the nuclei of galaxies. For this purpose, we formulate a power 

spectrum of the GWBR from cosmological evolving eccentric binaries. Then, we apply this 
formulation to the case of the GWBR from SMBH binaries. The key to doing this is to 
correctly estimate the number density of coalescing SMBH binaries. In this study, we use a 
semi-analytic model of galaxy and SMBH formation. We find that the power spectrum of 
the GWBR from SMBH binaries on eccentric orbits is suppressed for frequencies / < 1 nHz 
if the initial eccentricity, eo, satisfies eo > 0.2 and the initial semi-major axis is 300 times 
Scwarzschild radius. Our model predicts that while the overall shape and amplitude of the 
power spectrum depend strongly on the processes of galaxy formation, the eccentricity of 
binaries can affect the shape of the power spectrum for lower frequencies, i.e., / < 1 nHz. 
' Pulsar timing measurements, which can detect GW in this frequency range, could constrain 

' the effect of eccentricity on the power spectrum of the GWBR from SMBH binaries. 

o ■ 

^ ■ §1. Introduction 

O \ An ensemble of gravitational waves (GWs) from a number of inspiraling binaries 

-4— > ' of compact objects at different redshifts can be observed as stochastic gravitational 

\ wave background radiation (GWBR). Coalescing supermassive black hole (SMBH) 

binaries with masses in the rage 10^ - 10^ Mq emit GWs, results of the merging 
of their host galaxies. Recent studies^^'^-* predict that these SMBH binaries produce 
GWBR of frequencies ~ 1 n - 1 /iHz. In the frequency range 0.1 m - 10 mHz, it 
\ is believed that extragalactic close binaries of white dwarfs are dominant sources of 

GWBR. 3) 

Gravitational waves with frequencies in the range In- 100 nHz can be detected 
by pulsar timing measurements. Thus we should be able to detect the GWBR 
from SMBH binaries directly. Recently, Jenet et al.^^ developed a new method 
for detecting GWBR using multiple pulsar timing data. In the case of the Parkes 
Pulsar Timing Array (PPTA) project ,^^**^ their results showed that with regular 
timing observations of 20 pulsars with a timing accuracy of 100 ns, we should be 
able to directly detect the predicted levels of the GWBR from SMBH binaries within 
five years. ^) 

In previous studies^-*' of the GWBR from SMBH binaries, it was assumed that 
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all binaries are in circular orbits. However, binary orbits are generally eccentric. 
Before entering into the GW emitting regime, the evolution of a SMBH binary is de- 
termined by its dynamical interaction with field stars in the center of its host galaxy. 
Mikkola and Valtonen^) constructed approximate expressions relating the energy and 
angular momentum transfer rates for a heavy binary and the surrounding field of 
light stars, and they showed that the length of the semi-major axis decreases and the 
eccentricity increases owing to dynamical friction. Fukushige et al.^^ investigated 
the evolution of a SMBH binary in a uniformly distributed background of field stars 
and found that the dynamical friction on the eccentric binary is most effective at 
the apocenter, where the orbital velocity is minimum, and thus the eccentricity in- 
creases. Iwasawa ct al.^^ performed A^-body simulations of the dynamical evolution 
of triple SMBH (a binary of SMBHs and a single SMBH) systems in galactic nuclei 
and showed that the eccentricity of the SMBH binary reaches e > 0.9 in the GW 
emission regime. They found that both the Kozai mechanism and the thermalization 
of eccentricity due to the strong binary-single SMBH interaction drive the increase 
of the orbital eccentricity. Matsubayashi et al.^'^) investigated the orbital evolution 
of intermediate mass black hole (IMBH)-SMBH systems in galactic centers, using 
A?^-body simulations. They found that the eccentricity approaches unity (e > 0.8) 
and the IMBH can cause the SMBH to rapidly coalesce as a result of GW emission. 
Armitage and Natarajan^^) studied the evolution of the SMBH binary eccentricity 
which is excited by the interaction between a binary and circumbinary gas disk. They 
estimated a typical eccentricity at one week prior to coalescence to be e ~ 0.01. In 
the case of extreme-mass ratio binaries, higher eccentricities (e ~ 0.1) are possible. 

An eccentric binary emits GWs at all integer harmonics of the orbital fre- 
quency. -"^^^'^^^ Thus, the spectral energy distribution (SED) is different from that 
of a binary in a circular orbit, even if the masses and semi-major axes of the two 
binaries are the same. For larger eccentricities, radiation from higher harmonics are 
greater and thus the peak of the SED shifts toward higher frequency. Moreover, the 
orbital evolution of a binary due to GW radiation depends strongly on its eccen- 
tricity. As a result, the power of GW radiation of an eccentric binary is larger than 
that of a circular binary, and the timescale of GW radiation of an eccentric binary 
is shorter than that of a circular binary. Therefore, in order to predict the power 
spectrum of GWBR from compact binaries, it is necessary to take account of orbital 
eccentricities of binaries. For the cases of GWBRs from galactic and extragalactic 
neutron star binaries and black hole MACHO binaries, some studies have taken ac- 
count of the effect of eccentricity. ^^'^^^ However, for the case of the GWBR from 
SMBH binaries, none of models include the effect of eccentricity. 

In this study, we investigate the effect of orbital eccentricity on the expected 
GWBR from SMBH binaries. First, we formulate the power spectrum of GWBR 
from cosmological eccentric binaries, taking into account eccentricity evolution. In 
order to formulate the power spectrum, we adopt a simple relationship between the 
power spectrum of the GWBR produced by cosmological GW sources, the total 
time-integrated energy spectrum of individual sources, and the comoving number 
density of GW sources found by Phinney.^^^ Next, we apply this formulation to the 
case of the GWBR from coalescing SMBH binaries. In order to estimate the number 
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density of SMBH binaries, we use a semi-analytic (SA) model^''-* in which SMBH 
formation is incorporated into galaxy formation. 

This paper is organized as follows. In §2 we briefly review the GW emission of 
a binary in an eccentric orbit. In §3 we formulate the power spectrum of GWBR 
from cosmological compact binaries in eccentric orbits. In §4 we present the power 
spectrum of the GWBR from eccentric SMBH binaries. In §5 we present a summary 
and conclusions. 



§2. GWs from a binary in an eccentric orbit 



Here we briefly review the situation for GWs emitting from a binary in an 
eccentric orbit and the evolution of a binary due to GW emission in the weak field, 
slow motion limit. ^^^''^^^ 

The total power of the GW emission from a Keplerian binary consisting of two 
point masses Mi and M2, with orbital frequency fp and orbital eccentricity e is 

Lgw(Mi, M2, /p, e) = LGw,circ(Mi, M2, /p)F(e), (2-1) 

where 

32G^/3 



i^GW,circ(Mi,M2,/p) = :^^Mi°/J(27r/,)io/3 



and 



_ 1 + 73eV24 + 37eV96 
^(^) = (1 - e2)V/2 ■ (2-3) 

The quantity LGW,circ(^^i7 fp) is the total power from a binary in a circular orbit 
with masses Mi and M2 and orbital frequency fp. This is equal to the reciprocal 
of the proper rest-frame period of the binary. In the above expressions, G is the 
gravitational constant, c is the speed of light, and McMrp = [MiM2(Mi + M2)-^/^]^/^ 
is the chirp mass of the system. The distribution of the total power of GW emission 
among the harmonics of the orbital frequency is given by LGw,circ(-^i) -^2, fp)giji, e), 
with the rest-frame GW frequency fr = nfp. Here, g{n,e) is the GW frequency 
distribution function, expressed as 



4 

n 



g(n,e) = — 



2 



Jn-2{ne) - 2eJ„_i(ne) + -J„(ne) + 2eJ„+i(ne) - J„+2(rae) 



n 



+(1 - e') [Jn-2(ne) - 2eJ„(ne) + Jn+2ine)Y + ^ [Jn(ne)]"| , (2-4) 

where J„ is the nth-order Bessel function. It can be shown that 

00 

^g{n,e)=Fie). (2-5) 

n=l 
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For the case of a circular orbit (i.e., e = 0), wc have g{2, 0) = 1 and g{n, 0) = for 
all n 7^ 2. The SED of gravitational radiation is given by 



Lfrie, tp) = i^GW,circ(/p) XI Si^' ^)Kfr " "/p)- 



(2-6) 



n=l 



Here, tp is the time at which the orbital frequency is /p, and 5{x) is Dirac's delta 
function. 

The timescale of GW emission by a binary with orbital frequency fp measured 
in the rest frame is 

fv% (2-7) 



This timescale is given by 



TGw(Mi,M2,/p,e) = 



"^GW.circ 

F{e) 



(2-8) 



where TGw,circ(^i, Af2, /p) is the timescale of GW emission by a binary in a circular 
orbit. This timescale is given by 

t; / .3 \5/3 
96 VGMehirpJ 

■^rh i rn \ 2 fp ^ 



TGW,circ(Ml,-^2,/p) = — 



(27r/p)-«/3 



1.2 X 10^ 



o 



10^7 Hz 



yr. (2-9) 



As a result of the GW emission, the binary loses energy and angular momentum, 
and as a result, the orbit of the binary becomes more circular. The evolutions of the 
semi-major axis a and the eccentricity e are described by the follows, 

da _ _ 64 G^MiMsMtot 
'dt ~ 



73 2 37 4 

5 c5a3(l-e2)7/2 V ^ 24^ ^ 96^ 



64 G^MiM2Mtot 



de 
dt 



304 G^MiAhMt 



tot 



(2-10) 

15 c5a4(i _e2)5/2^^ > ^ ' ' ' ^^'-^^^ 

where Aftot = Mi + M2. Starting from a given orbit with parameters ao and cq, Eqs. 
(2-10) and (2-11) give the following relation between a and e: 



121 9 

■e 1 + 

304 



a 



1 



1 



12 

e ^ 19 
eo 



1 + i^e2 
^ ^ 304^ 

^ ^ 304^0 



(2-12) 



From the equality a'^ = G Mtot / {2tt fp)'^ , we find that the relation between the orbital 
frequency and the eccentricity is given by 

-3/2 



fp,0 



1 



1 



^ ^ 304 

121 „2 
304*^0 



1 + 



870 
2299 



(2-13) 
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where fp^ is the initial orbital frequency. In Fig. 1, we plot the evolution of the 
eccentricity as a function of the orbital frequency, /p//p,o- 




log[fp/fp,o] 

Fig. 1. Evolution of the eccentricity, e = e{fp/ fpfi,eo), as a function of the orbital frequency, 
fp/fp,o, for eo = 0.2, 0.4, 0.6, 0.8 and 0.9. 



§3. GWBR from binaries in eccentric orbits 

3.1. Characteristic amplitude of the GWBR spectrum 

In any homogeneous and isotropic universe, the present-day GWBR energy den- 
sity, PGwc?, must be equal to the sum of the energy densities radiated at each redshift 
z, divided by {1 + z) to account for the redshifting: 

PGWC^ = / / nc{z)—^^^^^dfrdz (3-1) 
Jo Jo i- + z ajr 

°°rn.(z)-^/.^d4. (3.2) 

7o 1 + 2^ dfr f 

Here, nc{z)dz is the comoving number density of GW sources at redshifts in the 
range z - z + dz. Also is the GW frequency in the source's rest frame, and / is the 
GW frequency in the observer frame^ and thus we have = f{l + z). The quantity 

dEQ\^ 

,r dfr (3-3) 

is the total energy emitted in GWs between the frequencies fr and + df,.. This 
energy is measured in the source's rest frame, and it is integrated over all solid 
angles and over the entire radiating lifetime of the source. Then, the total present- 
day energy density in GWs is 

PGWc'^l^^lp'hUlnf)^, (3-4) 
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where /ic(ln/) is the characteristic amphtude of the GWBR power spectrum over a 
logarithmic frequency interval dlnf = df/f. Therefore, we find 



hl(lnf) 



AG f°° , , 1 

T^C^f Jo 1 + 



AG 

TTC^f 



10 

roo 



+ z 



dEGw \ 

dfr J 



dz, 



roo 

Jo 



[dEGw] 
V dfr J 



fr=f{l + z) 

dz. 



fr = f{l+z) 



(3-5) 



If the GW sources are coalescing binaries, the characteristic amplitude is given by 



AG f 

hlilnf) = / dMidM2dz nc{MuM2,z) 



fdEGw{Mi,M2) 



dfr 



fr=f{l+z) 



.(3-6) 



where ndMi, M2, z)dMidM2dz is the comoving number density of coalescing bina- 
ries with masses Mi - Mi + dMi and M2 - M2 + dM2 (or z - z + dz. 

Note that Eq. (3-5) assumes that the timescale of GW emission is much shorter 
than the Hubble time. This relationship between the power spectrum of the GWBR 
produced by cosmological GW sources and the SED of an individual source was 
derived by Phinney.^^^ 

3.2. Total energy emitted in GWs 

The total energy of GW radiation emitted over a lifetime of duration tiHe is 

'•tliCe 



rtliie 

-E'gw = / LQSN{tp)dtp 
Jo 



Lfr{tp)dfrdtp, 



(3-7) 

(3-8) 



where Lgw(^p) is the power of GW emission and Lf^{tp) is the SED of GW emission. 
Then, we have 



dfr 



/•tlife 

/ Lf^{tp)dtp. 
Jo 



(3-9) 



Therefore, given the number density and SED of GW sources, we can calculate the 
amplitude of the GWBR power spectrum. 

3.3. Characteristic amplitude of the GWBR spectrum from eccentric binaries 

For a Keplerian binary consisting of two point masses Mi and M2 with an orbital 
eccentricity e, the SED of the GW radiation is given by Eq. (2-6). Prom Eqs. (2-6), 
(2-7) and (3-9), we obtain 



dfr 



I LGW,circ{fp)^ XI Si^' (^Wr " nfp)dfp 



= E 

n=l 



LGW,circ{fp 



nfp 



fp=fr/n 



(3-10) 
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Thus, from Eqs. (3-6) and (3-10), the power spectrum of GWBR from binaries in 
eccentric orbits is given by 



i'GW,circ(/p) — gyn, e) 



fp=f{l+z)/n 



AG f 

hl{\nf) = —— / dMidM2dznc{Mi,M2,z)Y^ 

J J n=l 

= dM,dM2dzn,{M,,M2,z){l + z)-y^ (^^^^ (^/)-4/3^, (3.II) 



where 

00 

^ = J2^n (3-12) 

n=l 

and 

The left panel of Fig. 2 plots g{n, e) as a function of the eccentricity, e, and the right 
panel plots <P and as functions of e. 

Here we note that the eccentricity is a function of the orbital frequency, fp/ fpfl, 
[see Eq. (2-13) and Fig. 1]. Therefore, e = e(/p//p,o, cq) = e{fr/nfpfi,eo) = 
e[/(l + z)/nfp^Q,eo]. Hence, (p depends on /, fp^, cq and z. Since SMBH bina- 
ries have various initial eccentricities, the number density in Eq. (3-11) should be 
nc{Mi, M2, eo, z) and it should be integrated over bq. If all binaries are circular (i.e., 
Co = 0), #2 = 1 and = for n 7^ 2. Thus, we have he oc /~^/^ for cq = 0. By 
contrast, for eg 7^ 0, # depends on /, and he is not proportional to /^^/^. Therefore, 
<P indicates the strength of the effect of orbital eccentricity on the GWBR power 
spectrum. In the left panel of Fig. 3, we plot <P as a function of rest-frame GW 
frequency, fr/fp,o, for some initial eccentricities. In the right panel of Fig. 3, we plot 
<P„ with Co = 0.8 as a function of fr/ fp,o for some harmonics. Owing to the radiation 
of the harmonics of fp, the power spectrum is suppressed at lower frequencies, i.e. 
fr/ fp,o ^ 10, and it is amplified at intermediate frequencies, i.e. 10 < fr/fp,o ^ 100. 
For larger frequencies, 100 < fr/fp,o, the main contributors to the power spectrum 
are circular binaries, which evolve from eccentric binaries and radiate only through 
the n = 2 mode. 



§4. Power spectrum of the GWBR from eccentric SMBH binaries 

4.1. A unified semi-analytic model of galaxy and SMBH formation 

In order to predict the power spectrum of the GWBR from coalescing SMBH 
binaries, we must estimate the number density of coalescing SMBH binaries. To do 
this, we use a semi-analytic (SA) model^^^ in which SMBH formation is incorporated 
into galaxy formation. 

In the standard hierarchical structure formation scenario in a cold dark matter 
(CDM) universe, dark-matter halos {dark halos) cluster through the effect of gravity 
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e e 

Fig. 2. Left panel: g{n,e) as a function of e for n = 1,2,3,4 and 5. Right panel: $ and for 
n = 1, 2, 3, 4 and 5 as functions of e. 




log[f/fp„] log[f/fp,o] 

Fig. 3. Strength of the effect of the orbital eccentricity on the GWBR power spectra. Left panel: 
$ = ^?(/r//p,o, eo) as a function of fr/fp,o for eo = 0.2,0.4,0.6,0.8 and 0.9. Right panel: 
with eo = 0.8 as a function of /r//p,o for n = 1, 2, 3, 4 and 5. $ for eo = 0.8 is also shown. 

and merge. In each merged dark halo, a galaxy is formed as a result of radiative 
gas cooling. In each galaxy, star formation and supernova feedback occur. Several 
galaxies in a common dark halo sometimes merge, and a more massive galaxy is 
thereby formed. When galaxies merge, SMBHs in the centers of these galaxies move 
toward the center of the new merged galaxy and subsequently form a SMBH binary. 
If the binary loses a sufficient amount of energy and angular momentum, it will evolve 
into the GW emitting regime and begin inspiraling, eventually coalescing with a GW 
burst. 

In SA models, merging histories of dark halos are realized using a Monte-Carlo 
algorithm, and the evolution of baryonic components within dark halos is calculated 
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using simple analytic models for gas cooling, star formation, supernova feedback, 
galaxy merging and other processes. SA models have successfully reproduced a 
variety of observed features of galaxies, such as their luminosity functions and gas 
fractions in disk galaxies. In our SA model in which SMBH formation has been 
incorporated,^^) it is assumed that SMBHs grow through the coalescence of two or 
more SMBHs when their host galaxies merge. We also assume that during a major 
merger, a fraction of the cold gas proportional to the total mass of stars newly formed 
at the starburst is accreted onto the newly formed SMBH and that this gas fueling 
leads to quasar activity. Thus, SMBHs also grow through the accretion of cold gas. 
Our SA model reproduces not only observational features of galaxies^^) but also the 
present-day observed SMBH mass function and the quasar luminosity functions at 
different rcdshifts.-'^'^) Using this SA model, Enoki et al. estimated the amplitude of 
GWBR from inspiraling SMBH binaries in circular orbits and the event rate of GW 
bursts due to SMBH binary coalescing.^) 

It is difficult to determine how SMBH binaries manage to shrink into a regime 
in which their evolution is driven by GW emission after their host galaxies merge, 
because not all the physical processes and conditions related to this problem (dy- 
namical friction, the stellar distribution, triplet SMBH interaction, gas dynamical 
effects, and so on) are yet clear (see, e.g. Refs.'^)'^) and"^^)). Therefore, it is also 
difficult to predict the distribution of the initial eccentricity. In what follows, for 
simplicity, we assume that all SMBH binaries coalesce when host galaxies merge. In 
other words, the efficiency of SMBH coalescence is assumed to be maximal. Thus, 
the efficiency of SMBH coalescence is maximal and the predicted amplitude of the 
GWBR spectrum should be interpreted as the upper limit (see Enoki et al.^-*). In 
this paper, in order to clarify the effect of eccentricity, we assume that all binaries 
have a single initial eccentricity given at a radius that is equal to some constant 
times the Schwarzschild radius of the SMBH. 

In this study, the adopted cosmological model is a low-density, spatially flat 
cold dark matter (ylCDM) universe with the present density parameter, = 0.3, 
the cosmological constant J7yi = 0.7, and the Hubble constant h = 0.7 (h = 
Hq/100 km s^^ Mpc^^). A detailed description of the model and its parameters 
is given in Nagashima et al.^^-* and Enoki et al."'^)'"^''') 

4.2. Maximum orbital frequency 

Substituting the number density of coalescing SMBH binaries, nc{Mi, M2, z)dMidM2dz, 
into Eq. (3-11), we can calculate the power spectrum of the GWBR from SMBH 
coalescing SMBH binaries. Here, we introduce the cut off orbital frequency /p,max-^^ 
As a binary evolves with time owing to its emission of GW, the frequency increases. 
We assume that the binary orbit is quasi-stationary until the radius equals 3i?s, 
where Rs is the Schwarzschild radius, i.e. the radius of the innermost stable circular 
orbit (ISCO) for a particle around a non-rotating black hole. Then, the maximum 
orbital frequency /p,max is 

c^ / MoX^/^ 
WM,M,) = j^3^j^-^(l + -j 
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= 2.2 X 10-^ 1 + — Hz, (4-1) 

where Mi and M2 are the SMBH masses (and we assume Mi > M2). Then, in Eq. 
(2-6), we replace -LGW,circ(/p) by LGW,circ(/p)6'(/p,max - /p), where e{x) is the step 
function. Therefore, becomes 



= (^) (/P,max - /(I + ^)/n) . (4-2) 



If this cut off does not exist, we have he (x f for / > 2/p^inax- 
4.3. Results 

In the left panel of Fig. 4, we plot power spectra of GWBR, he, from SMBH 
binaries for various initial eccentricities. The initial eccentricities, eg, are those 
at /p,o//p,max = 10^^ (a = 300i?s)- This figure shows that power spectra at the 
frequencies measured by pulsar timing (~ 1 n — 100 nHz) are suppressed as a result 
of harmonic radiation, especially for cq ^ 0.4. The slope of the spectrum changes 
near / = 1 /^Hz, owing to a lack of power associated with the upper limit frequency, 

/p,max- 

The right panel of Fig. 4 displays power spectra of GWBR from SMBH binaries 
for Co = 0.8 for different initial orbital frequencies with /p,o//p,max = 10~^, 20~^ 

_q —2/3 

and 30~ . Because a oc fp , these correspond to a = 75 i?s,300 Rs, 1200 i?s and 
2700 Rs, respectively. Even in the case /p,o//p,max = 30~^, the power spectrum for 
/ < 1 nHz is suppressed. We note that for a binary with cq = 0.8 at /p,o//p,max = 
30^^ (a = 2700i?s), the eccentricity of the binary is e ~ 0.95 for a ~ 10^ Rs ~ 
0.1 (Mbh/IO^Mq) pc. Both panels of Fig. 4 indicate that the shape of the power 
spectrum in the low frequency range depends strongly on eo and /p,o//p,max- 

In Fig. 5, wc plot some harmonics of the GWBR power spectrum from SMBH 
binaries for cq = 0.8 and /p,o//p,max = 10~^. As with the right panel of Fig. 3, 
in the lower frequency range, / ^ 1 nHz, the n = 1 mode is dominant and in the 
higher frequency range, / ^ 0.1 fjHz, the n = 2 mode is dominant. Therefore, in 
this frequency range, GWs from circular binaries are dominant. 

The present-day energy density parameter of GWBR over a logarithmic fre- 
quency interval d In / = d/// is given by 

f2Gw(lnf) = ^fhl{lnf). (4-3) 

In Fig. 6, we plot /?Qw^(lny) from SMBH binaries with 6q — 0.8 and /p.o/ /p.max — 
10~^ and with eo = for several intervals of the total SMBHs mass (Mtot = M1-I-M2). 
It is seen that the main contribution to the energy density is from binaries with total 
masses Mtot = 10^ — 10^*^ Mq. In the lower frequency range (/ < 0.1 /xHz), which 
can be measured by pulsar timing, GWs from binaries with masses Mtot > lO^M© 
are dominant. As with the left panel of Fig. 3, Fig. 6 shows suppressions and ampli- 
fications of the energy density parameters for each SMBH mass interval due to the 
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log(f[Hz]) log(f[Hz]) 

Fig. 4. Characteristic amplitudes of GWBR power spectra over a logarithmic frequency interval, 
hc{lnf), from SMBH binaries. Left panel: Power spectra of GWBR from SMBH binaries with 
/p,o//p,max ~ 10~^ for several initial eccentricities, eo = 0.0,0.2,0.4,0.6,0.8 and 0.9. Right 
panel: Power spectra of GWBR from SMBH binaries for eo = 0.8 for several initial orbital 
frequencies, ,fp,o/,fp,max = 5"'^, 10"'^, 20"'^ and 30"'^. In each panel, the solid triangle indicates 
the current limit from pulsar timing measurements.^^-* 




-10 -9 -8 -7 -6 



log(f[Hz]) 

Fig. 5. Power spectra of GWBR from SMBH binaries for eo = 0.8. Some harmonics of the GWBR 
spectrum for n = 1,2 and 3. The initial orbital frequency is /p,o//p,max = 10""^. We also plot 
the power spectrum of GWBR from circular binaries (dashed line) and the current limit from 
pulsar timing measurements (solid triangle). '^^^ 

harmonic radiation. Here, the initial eccentricity is that at /p.o//p.max = 10~^. From 
Eq. (4-1), we find fp^ cx /p,max oc (MiMtot)"^''^- Thus, for smaller-mass SMBHs, 
non-zero eccentricities affect the power of higher frequencies. When we sum up the 
energy densities of GWBR from SMBH mass intervals at approximately 10 nHz, 
where the energy density of GWBR from SMBH binaries with Mtot > 10^ Mq in- 
creases, the energy density of GWBR from SMBH binaries with Mtot < 10^ Mq has 
decreased significantly. For this reason, we cannot observe the amplification of the 



12 



M. Enoki and M. Nagashima 



total energy density of GWBR from SMBH binaries. 

Pulsar timing measurements allow GWs with frequencies in the range 1 n — 
100 nHz to be detected. Under the assumption that he oc Z^^^^, the completed 
PPTA data set (twenty pulsars with an rms timing residual of 100 ns over 5 years) 
could potentially provide limits on the level of GWBR from SMBH binaries: f^cw = 
5.5 X 10-1° at / = 10-^-5 Hz (1 yr'^), Oqw = 1.3x10-^° at / = IQ-^-^ Hz (1/8 yr'^) 
and new = 7.3 x lO'^^ at / = lO"^-^ Hz (1/20 yr-i).20) In Fig. 6, we also plot the 
potential future limits. These limits are well under the predicted energy density of 
GWBR from SMBH binaries for eo = {he oc /"^Z^). Therefore, we conclude that 
the full PPTA data set should be sufficient to constrain the effect of eccentricity on 
the GWBR from SMBH binaries. 
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Fig. 6. Energy density parameters of GWBR over a logarithmic frequency interval f2Gw(ln/) in 
several intervals of the total mass. The thick curves represent the results for eo — 0.8 and 
/p,o//p,max ~ W~^, and thin curves represent the results for eo — 0. The solid triangle indicates 
the current limit from pulsar timing measurements.'^^-' The three other triangles indicate the 
potential future lower limits from the full PPTA data set for the case he oc j-2/3 20) 



4.4. Effects of galaxy formation processes on GWBR from SMBH binaries 

The number density of coalescing SMBH binaries with masses in the range Mi — 
Ml + dMi and M2-M2 + dM2 at z-z + dz, i.e. n^Mi , M2 , z)dMi dM2dz, depends 
on galaxy formation processes. In our SA model, the dominant mass growth process 
of SMBHs is the accretion of cold gas, which is also the material of stars. Thus, 
nc{Mi, M2, z) depends strongly on processes related to star formation related. Here, 
we show how the star formation time scale and the strength of supernova feedback 
affect the power spectrum of the GWBR from SMBH binaries. 

In our SA model, the star formation rate (SFR) of a galaxy is assumed to 
be M* = Mcom/t*, where Mcqm is the mass of the cold gas and r* is the time 
scale of star formation. We assume r^, = T^(l/circ/300 kms"^)"*, where Vdrc is 
the circular velocity of the galaxy. The free parameters and a* are chosen to 
match the observed mass fraction of cold gas in the disks of spiral galaxies. With 
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star formation, supernovae occur and heat up the surrounding cold gas, yielding a 
hot gas phase (a process called supernova feedback). The reheating rate is given 
by Mrcheat = /?(V'circ)M*, whcrc [3{Vc\tc) = (Vhot /V'circ)"''"' • The free parameters 
T4iot a-iid ajiot are determined by matching the observed local luminosity function 
of galaxies. In a previous study^^ (and in a previous subsection of this paper), we 
adopted t° = 1.5 Gyr, a* = —2, l^iot = 280 km s~^ and Ohot = 2.5 as fiducial values. 

The left panel of Fig. 7 displays power spectra of GWBR from SMBH binaries for 
various star formation time scales, r^. Other parameters are the same in each case. It 
is seen that for large r^, the SFR is small and thus a large amount of cold gas remains 
in each galaxy. Thus, the masses of SMBHs become large. Therefore, the amplitude 
of the power spectrum of the GWBR becomes large. Moreover, the slope changing 
frequency moves forward lower frequency, because the upper limit frequency, /p^maxj 
becomes small. The right panel of Fig. 7 plots power spectra of GWBR from SMBH 
binaries for various supernovae feedback strengths, Vhot- Other parameters arc the 
same in each case. In the case of no supernovae feedback (I^ot = km s^^), the 
cold gas is not heated to a hot gas. Thus, a large amount of cold gas remains in 
each galaxy, and the masses of SMBHs become large. By contrast, for the large Vhot 
case, the cold gas is greatly heated, and thus the masses of SMBH become small. 
Both panels of Fig. 7 indicate that although the shape of the power spectrum of 
the GWBR from SMBH binaries is affected by the eccentricity of binaries for low 
frequencies, / < 1 uHz, the shape of the power spectrum for higher frequencies and 
the overall amplitude depend strongly on the processes of galaxy formation. 

We note that the purpose of this subsection is to demonstrate the importance 
of galaxy formation processes on the growth of SMBHs and the GWBR from SMBH 
binaries. Therefore, in this subsection, we present cases with extreme parameter 
values. Results for models with these parameter values are greatly inconsistent with 
observational results, such as galaxy luminosity functions and the present SMBH 
mass function. In Fig. 8, we plot SMBH mass functions at z = 0. It is seen 
that only the model with fiducial values is consistent with the observational results 
obtained by Salucci et al.^^^ 

§5. Summary and conclusions 

In this study, we have investigated how orbital eccentricities of binaries affect 
power spectra of GWBR from coalescing SMBH binaries. A compact binary in 
an eccentric orbit radiates GWs at all integer harmonics of its orbital frequency. 
Owing to this harmonic radiation, the SED, the power and the timescale of the GW 
emission of a binary in an eccentric orbit are different from those of a binary in a 
circular orbit. Therefore, first, we formulated the power spectrum of GWBR from 
cosmological compact binaries in eccentric orbits. Then using this formulation and 
our SA model for galaxy and SMBH formation, we calculated the power spectra of 
GWBR from coalescing SMBH binaries in eccentric orbits. 

We found that the calculated power spectra of the GWBR from SMBH binaries 
in eccentric orbits are suppressed owing to the harmonic radiation for lower frequen- 
cies (/ < 1 nHz) if the initial eccentricity satisfies eo > 0.2 at a = 300 i?s- The 
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Fig. 7. Effects of galaxy formation processes on the power spectrum of GWBR from SMBH binaries 
with eo = 0.8 and /p,o//p,max = 10~^. Left panel: For three star formation time scales, 
Tt.o = 0.5, 1.5 and 4.5 Gyr. Here r,.o = 1.5 Gyr is the fiducial value. Right panel: For three 
supernovae feedback strengths, Vhot = 0,280 and 420 km s^^. Here, Vhot = 280 km s^^ is the 
fiducial value. The case Vhot = km s~^ corresponds to no supernova feedback. In each panel, 
we also plot the current limit obtained from pulsar timing measurements (solid triangle). 




log[MBH/MsuN ] 



Fig. 8. SMBH mass functions of several models at 2: = 0. The solid curve represents the results of 
the model with fiducial values, t*,o = 1-5 Gyr and Vhot = 280 km s~^. The dots with error bars 
present the data for the present mass function obtained by Salucci et al.^^^ 



degree of this suppression depends strongly on the initial eccentricity. In this paper, 
for simplicity and in order to clarify the effect of the eccentricity, we assumed that 
all binaries have the same initial eccentricity. The results of our study suggest that 
elucidating the initial eccentricity distribution is essential for the investigation of the 
power spectrum of GWBR from compact binaries. 

Because the number density of coalescing SMBH binaries is determined by the 
processes of galaxy formation, the power spectrum of the GWBR from SMBH bina- 
ries depends strongly on these processes, especially processes related to star forma- 
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tion, which regulate the amount of cold gas accreted onto SMBHs. The prediction 
of our model shows that the overall shape and amplitude of the power spectrum of 
the GWBR from coalescing SMBH binaries depend on galaxy formation processes. 
However, for low frequencies (/ ^ 1 nHz), the shape of the power spectrum of the 
GWBR also depends on the initial eccentricity. Because the pulsar timing measure- 
ments can provide limits on the level of GWBR with frequencies / ~ 1 n — 100 nHz, 
pulsar timing observations, such as PPTA project, should provide constraints not 
only on the number density of coalescing SMBH binaries but also on the effect of 
orbital eccentricity on the GWBR from SMBH binaries. 
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